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ABSTRACT 

In this paper we analyze the behavior of Galactic, LMC and SMC Cepheids in terms 
of period-color (PC) and amplitude-color (AC) diagrams at the phases of maximum, 
mean and minimum light. We find very different behavior between Galactic and Mag- 
ellanic Cloud Cepheids. Motivated by the recent report of a break in LMC PC relations 
at 10 days (Tammann et al. 2002), we use the F-statistical test to examine the PC 
relations at mean light in these three galaxies. The results of the F-test support the ex- 
istence of the a break in the LMC PC (mean) relation, but not in the Galactic or SMC 
PC(mean) relations. Furthermore, the LMC Cepheids also show a break at minimum 
light, which is not seen in the Galactic and SMC Cepheids. We further discuss the 
effect on the period-luminosity relations in the LMC due to the break in the PC (mean) 
relation. 
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1 INTRODUCTION 



ISimon. Kanbur fc MihalasI (Il99l hereafter SKM) used 
hydro -dynamical models to explain the observations of lCodd 
l|l943): Galactic Cepheids show a spectral type independent 
of period at maximum light and a spectral type at min- 
imum light that gets later as the period increases. SKM 
computed radiative hydro-dynamical models of Galactic 
Cepheids which agreed with Code's observations. SKM in- 
terpreted these observational phenomena as being due to 
the location of the photosphere relative to the hydrogen ion- 
ization front. They further used the Stefan-Boltzmann law 
appUed at maximum and minimum light to show that 



logT™ 



\ogTrain = YO^^™ 



(1) 



where Tmax and Tmin refer to the effective temperature 
at maximum and minimum light, respectively. Thus, higher 
optical amplitudes are associated with higher temperature 
amplitudes, which are in turn related to higher values of 
Tmax and /or Tmin- For this study we do not assert a causal 
relation between higher temperatures and higher amplitudes 
preferring to refer to these quantities as being "associated" . 
If, for some reason, either Tmax or Tmin does not increase 
as the amplitude increases, equation (1) predicts a relation- 
ship between am plitu d e and Tmin or Tm ax resp ectively . SKM 
used data from |Pe3 lll97d ) and iMoffctt fc Barnes! l|l98d 
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Il984f) to show that Galactic Cepheids are such that higher 
amplitude stars are driven to cooler temperatures at mini- 
mum light. This, according to equation (1), is because the 
range of temperatures at maximum light is independent of 
period for a large range of periods. So the form of the period- 
color (PC) relation at maximum light is related to the form 
of the amplitude-color (AC) relation at minimum light, and 

vice versa^ 

Tam mann et al.l ll2003ll used Galactic and OGLE 
LMC/SMC Cepheid data to show that there is a dif- 
ference in the PC relations in t hese three galaxies. Fur- 
thermore, iTammann et al.l (l2002l) and iTammann fc ReindJ 
li2002lh show the existence of two PC relations, one for short 
{P < 10 days) and one for long (P > 10 days) period 
Cepheids, in the LMC. Motivated by this and the presence 
of high quality Magellanic Cloud C epheid data from the 
OGLE project iUdalski et alJ[l999a^ ■ we decided to investi- 
gate the properties of Magellanic Cloud Cepheids in terms 
of their PC and AC relations at maximum, mean and mini- 
mum light. 

In addition to the two major reasons mentioned above, 
we list a number of other arguments motivating the present 
study: 

(i) Since mean light is just that - the average over a range 
of values - interesting properties of Cepheids at mean light 
are due to the average of these properties at all pulsation 
phases. By investigating the phases of maximum and mini- 
mum light, we are studying those phases of stellar pulsation 
which contribute to the observed properties at mean light. 
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Our interest lies in understanding breaks in the LMC mean 
light Cepheid PC and period-luminosity (PL) relations at 10 
days reported bv lTammann et "aJ] (|22o2)- Let yij be the (ab- 
solute) magnitude of Cepheid variable stars, i = 1, ...,nstar 
at the j'*'' phase (j — 1, A'^) during a pulsation period Pi. 
Then we can formulate a PL relation at a particular phase 
as, 

yij ^ aj +hj\og{Pi), (2) 

where Uj , hj are the unknown coefficients as a function 
of the phase j. If we define yi as X/j-i Vij/^y ^^e av- 
erage over the pulsation period, it is easy to show that 
yi = a.m + bmlog{Pi), where am,bm are the average over 
phase of the aj,bj in equation (2)^. This will be true if y 
is measured in intensities and then the intensity mean con- 
verted to magnitudes or if y is measured in magnitudes. 
Of course the magnitude mean and intensity means are in 
general not equal to each other but the difference is small 
(~ O.OSmag. ) and constant ove r a wide range of periods (for 
example, see lGieren et al .Il998l) . Consequently, some insight 
into the mean light relation can be gained by studying PL 
relations at various phases, e.g. at maximum and minimum 
light. Since the PC relation affects the PL relation (see, e.g., 
[Madorc & Frcodman 1991 for the basic physics of PC and 
PL relations), it is of interest to study the PC relation at 
various phases. Furthermore, the maximum and minimum 
light are closely associated with the more interesting phases 
of stellar pulsation: the expansion/contraction velocity is 
close to its maximum value when the photosp here is passing 
through the mean radius of the star (see, e.g.. iMihala jl2003 
for the details). 

(ii) The amplitude is a fundamental observational and 
theoretical quantity in stellar pulsation. Kanbur and Ngeow 
(2004, in-preparation) show that the amplitude is a very 
good descriptor of light curve shape, and the V band ampli- 
tudes are correlated to the first Principal C omponent (PCI) 
of the light curve. In addition. iKanbur et al. (2002) demon- 
strated that PCI can explain over 90% of the variation in 
light curve shape. Thus the V band amplitude can be taken 
to be a good descriptor of V band light curve shape. A sim- 
ilar conclusion holds for the I band. Because the optical 
brightness fluctuations of Ce pheids are predominantly due 
to temperature fluctuations l|coxlll98nh . it is thus instruc- 
tive to examine AC diagrams at maximum and minimum V 
band light. Furthermore, since AC relations are related to 
PC relations through equation (1), their study can serve as 
a useful complement to strengthen any conclusions reached 
using PC relations. 



2 THE DATA 

The standard Johnson- Cousins V and I band data for the 
fundamental mode Cepheids in the Galaxy, LMC and SMC 
were used in this study. For both of the LMC and SMC 
Cepheids, the periods, V and I band photometric data 

^ Note that bm may not lie in between bmax and femin > the slopes 
at maximum and minimum light, respectively. A similar conclu- 
sion also holds for the zero-point, am ■ 



and the E{B — V) values for every Cepheid were taken 
from the OGLE JUdalski et al.ll999bll3l website^. There are 
771 and 488 fundamental mode Cepheids (as classifled by 
OGLE team) in LMC and SMC, respectively. For Galactic 
Cepheid data, we considered only the Cepheids classified as 
"DCEP" in t he General Cata log of Variable Stars (GOVS, 
fourth edition, lKholoDovll993) . The periods for the Galactic 
Cepheids are taken from the McMaster Cepheid dat abase^, 
and th e E{B — V) values are adopted from Ta mmann et alJ 
J2003l listed as E{B - V)corr in their table 1). The pho- 
tometric data for Gala ctic Cepheids w ere obtained from 
two sources: (a) Moffctt fc Earned ill 984^ . where actual data 
we re downloaded from the McMaster Cepheid da tabase; and 
(h) iBerdnikwl lll997t) '* and lBerdnikoy &: Turned i200n . For 
Cepheids that have e ntries in both iBer^fikovf Jl997li and 
iBerdnikov fc Turner! ||200 j) . photometric data from these 
two sources were merged together to provide a better light 
curve. Since the bandpasses used in iMoffett fc Earned (ll984^ 
are in Johnson V and I, we converted the Johnson I band 
photometric data in this datas et to Cousins I band with the 
color transformations given in ICoulson et al] ^ll985^ . 

The Cepheid photometric data in these three galaxies 
were then fitted with a Fourier expansion of the followin g 
form JSchaltenbrand fc Tammannll97ll:lNgeow et al.l2003f) : 

M 

m{t) = mo + ^[A^cos{2ni^{t) + (l>i)], (3) 

i = l 

where $(t) = {t — to)/P~'int[{t~to)/P], with to being a 
common starting epoch for all Cepheids in both bands. We 
used a simulated annealing technique to fit the da t a wit h 
this Fourier expansion, as described in lNgeow et alJ i2003l) . 
Thus our Fourier fits are carried out to published photomet- 
ric data and our means are magnitude means. For Galac- 
tic data, we adopted a fifth order expansion {M = 5) to 
most of the Cepheids. However, in certain cases a fourth or 
sixth order Fourier expansion gave a better fit to the data. 
For OGLE LMC Cepheids, w e fit the data with M = 4, 
as this dataset is a lso used in iNgeow et alJ (l2003h and in 
iKanbur et all jioQi). For OGLE SMC Cepheids, we mainly 
fit fourth or fifth order Fourier expansions to the data, while 
some of them were fitted with sixth order. All fitted light 
curves were then visually inspected (se e Kanbur et a l . 2003 
for the selection criteria). Figure 7 of iNgeow et al.l (|20o3) 
illustrates the improvement that can be obtained when us- 
ing our fitting method to the observed light curves. We 
only selected those Cepheids with well fitted light curves 
in both V and I bands. In addition, we exclude Cepheids 
with log(P) < 0.4 to avoid possible conta mination from flrst 
overtone Cepheids dUdalski et aLlll999al). The numbers of 
Cephe ids in the flnal samples are: 79 from lMoffett fc Earned 
lia&A data; 75 from Berdnikov data; 634 from OGLE LMC 
data; and 391 from OGLE SMC data. From the Fourier flts 
we obtained the following quantities: 

(i) V and I band amplitude from the numerical maxi- 
^ http : //bulge . pr inceton . edu/'^ogle/ 

^ http : //dogwood . physics . mcmaster . ca/Cepheld//HomePage . html 

^ Note that the column oi R — Ic should be swapped with column 
oi V — Ic in this dataset. 
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mum and minimum of the Fourier fits: Vamp = Vmin — Vmax, 

lamp — Imin ^max- 

(ii) (V^ — I)max: defined as Vmax — Iphmax, where Iphmax 
is the I band magnitude at the same phase as Vmax- 

(iii) {V — I)mean- defined as Vmean - Iphmean, where 
Iphmean IS the I baud magnitude at the same phase as Vmean- 
Vmean IS the V band magnitude closest to mo, the mean 
value in equation (3) . 

(iv) {V - I)min- defined as Vmin - Iphmin, where Iphmin 

is the I band magnitude at the same phase as Vmin- 

Finally, the colors at these three phases have been cor- 
rected for extinction using Av,i = Rv,iE{B — V). The val- 
ues of R are: Ry = 3.f7, Ri = f.89 for Galactic data 
jTammann et al.1 12003'). a nd Ry = 3.24 Ri = f.96 for 
OGLE LMC and SMC data iUdalski et al.l lll999blcl) . In fact, 
the results are unchanged for the values o f Ry and Ri as 
long a s AT? = Ry — Ri = 1.28 as given in lTammann et all 
l|2003h . because (y~I)o = {V~I)~ARE{B-V)- We apply 
the same extinction values of Ayj to the colors at maximum, 
mean and minimum light, since the quantities of Ayj should 
remain unchanged for any pulsation phases. In addition, our 
results are unchanged if we adopt {< V > — < I >)o for 
the mean color as in item (iii) above. 



3 ANALYSIS AND RESULTS 

Figures0|5|and|3|show the plots of log(P) vs. {V ~I)o color 
and Vamp vs. {V — I)o color at maximum, mean and mini- 
mum light (see the definition in Section 2) for the Cepheids 
in Galaxy, LMC and SMC respectively. Open and closed 
circles in these figures are for short (log(P) < 1.0) and long 
(log(P) > 1.0) period Cepheids, respectively. In what fol- 
lows short and long period Cepheids are always used in this 
sense. For all three data sets in the Galaxy, LMC and SMC, 
we can fit linear relations between log(P) and {V ~I)o color 
and Vamp and {V — I)o color, both for the entire sample, 
and also for short and long period Cepheids separately. Ta- 
bles and 121 show these results for the fits to period-color 
and amplitude-color relations, respectively. In these tables, 
column 1 gives the phase at which the fit is made - either 
at maximum, mean or minimum light. Column 2 is labelled 
MRSS and represents the mean residual sum of squares 
from the fit to the entire sample. The slope of the linear 
fit together with its associated error is given in column 3. 
Columns 4 and 5 give the residual mean sum of squares and 
slope plus error for the short period Cepheids. Columns 6 
and 7 do this for the longer period Cepheids. We discuss 
columns 8-11 shortly. In Figures fH3l we have plotted the fit- 
ted lines for short and long period Cepheids as solid and 
dashed lines respectively. 

What is of interest for the present study is whether these 
PC plots show statistical evidence of a change of slope at a 
period of 10 days. To investigate this, we can fit a regression 
line over the entire period range and then fit two regression 
lines, one separately for short and long period Cepheids. The 
former case is the reduced model while the latter case with 
two separate regression lines is the full model. We can write 
the model under consideration as 



Here Y is the dependent variable, in our case {V — I)o 
color at any of the three phases under consideration. Ws is 
an indicator variable which is 1 if the star's period is less 
than 10 days and otherwise. Wl is an indicator variable 
which is if the star's period is less than 10 days and 1 oth- 
erwise. The variable Zs contains the independent variable, 
either log(P) or Vamp, but is zero if the period is greater 
than 10 days. Zl is similar but is zero for periods less than 
10 days. The parameters as, aL, bs and &l are the zero- 
points and slopes for short and long period Cepheids, respec- 
tively. Thus what is of interest is if the data are consistent 
with bL = bs'- the s l ope is same for long and short period 
Cepheids. IWeisber'3 il980h shows that in this situation an 
appropriate F test statistic can be formulated as described 
in the following equation. 



F = 



{RSSr - RSSF)/[{n - 2) - (n - 4)] 
RSSp/in-A) 



(5) 



where RSSr,RSSf are the residual sums of squares 
in the reduced (single line regression) and full model (two 
lines regression), respectively, and n is the number of stars 
in the entire sample. We refer the F statistic in equation 
(5) to an F distribution, F2,n-4, under the null hypothe- 
sis that the two-parameter regression (i.e. a single line) is 
sufficient. The four-parameter regression (i.e. two lines) will 
have a smaller residual sum of squares. The probability of 
the observed value of the F statistic, P{F), under the null 
hypothesis, gives the significance of this reduction in sum of 
squares. Thus, a "large" value of F indicates that the null 
hypothesis can be rejected. Column 8 of TablesQand|21gives 
the values of the F statistic and column 9 gives its proba- 
bility value as described above. Note that because there are 
fewer observed long period Cepheids, the error on the slope 
for the long period PC relations is generally larger but this 
is automatically taken account of by the F test. 

It could also be that the short and long period slopes 
are similar but that the F statistic produces a significant 
result because the intercepts are different. In this case, we 
can also compare the statistical significance of differences in 
the slopes of these regressions by referring the quantity. 



t = {bs 



Sbr 



MS 



MS _^ MS 
SSs SSl ' 



SSs + SSl 



ns + riL 



(6) 

(7) 
(8) 



Y = asWs + ffliWi + bsZs + + e. 



(4) 



to a t distribution with ns -\- tll — ^ degrees of freedom. 
Pit) is the probability of the observed value of the t statistic 
under the null hypothesis that the two slopes are equal. In 
the above formulae SSs, SSL,ns,nL,bs,bL are the residual 
sum of squares, the number of Cepheids and the slopes in 
the PC or AC relation for short and long period Cepheids, 
respectively. Columns 10 and 11 in Tables and |5| give the 
values of t statistic and its associated P{t) value from two- 
tail t distribution. 

Our figures do not contain error bars on the estimated 
values of the {V — I)o color. These error bars are typically 
less than the size of the symbol. This includes errors due 
to photometry and reddening. A typical V band photomet- 
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Figure 1. The period-color and amplitude-color relations at phases of maximum, mean and minimum light for Galactic Cepheids. The 
open circles and filled circles correspond to short and long period Cepheids, respectively. The dashed and solid lines are the fitted period- 
color and amplitude-color relations to short and long period Cepheids, respectively. The number of short and long period Cepheids are 
^short = 113 and nia„g = 41. Left (a): The period-color relations. Right (b): The amplitude-color relations. 
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Figure 2. The period-color and amplitude-color relations at phases of maximum, mean and minimum light for OGLE LMC Cepheids. 
The open circles and filled circles correspond to short and long period Cepheids, respectively. The dashed and solid lines are the fitted 
period-color and amplitude-color relations to short and long period Cepheids, respectively. The number of short and long period Cepheids 
are nghart = 585 and niang = 49. Left (a): The period-color relations. Right (b): The amplitude-color relations. 



ric error quoted by the OGLE team is 0.015rnags and per- 
haps as much as 0.05mags for reddening. If we add these in 
quadrature, a rough estimate of sigma for the photometry is 
0.052mags. If we assume a similar figure for the I band (with 

errors of ~ O.OOlmags. to 0.005mags) and again add up the 
errors in quadrature, then a typical error for the {V — I)o 



color is ~ 0.07mags. Strictly speaking this error should be 
included in a weighted least squares fit when constructing 
our period-color or amplitude-color fits. The weights are the 
inverse of the error assigned to each data point. However if 

the weights arc constant for each star then they cancel out 
in the calculation of the F statistic and in the derivation of 
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Figure 3. The period-color and amplitude-color relations at phases of maximum, mean and minimum light for OGLE SMC Cepheids. 
The open circles and filled circles correspond to short and long period Cepheids, respectively. The dashed and solid lines are the fitted 
period-color and amplitude-color relations to short and long period Cepheids, respectively. The number of short and long period Cepheids 
are ngi^^rt = 334 and ni^ng = 57. Left (a): The period-color relations. Right (b): The amplitude-color relations. 



Table 1. Period-color relations at maximum, mean and minimum light. 



Phase 






MRSS'i; 




MRSSl 


Slopel 


F 


P(F) 


t 


P(t) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


Galactic 


Max 


0.009 


0.142 ±0.031 


0.006 


0.389 ± 0.062 


0.013 


0.020 ±0.105 


7.95 


0.00 


3.40 


0.00 


Mean 


0.008 


0.264 ±0.030 


0.006 


0.229 ± 0.062 


0.013 


0.373 ±0.104 


1.05 


0.35 


1.34 


0.18 


Min 


0.007 


0.310 ±0.028 


0.005 


0.330 ±0.058 


0.013 


0.309 ±0.104 


0.06 


0.94 


0.21 


0.83 


OGLE LMC 


Max 


0.010 


0.205 ±0.020 


0.010 


0.234 ± 0.030 


0.008 


-0.031 ± 0.101 


2.67 


0.07 


2.24 


0.03 


Mean 


0.009 


0.239 ±0.018 


0.008 


0.152 ±0.027 


0.007 


0.590 ± 0.097 


13.19 


0.00 


4.10 


0.00 


Min 


0.008 


0.291 ±0.018 


0.008 


0.204 ±0.027 


0.009 


0.493 ±0.109 


10.38 


0.00 


2.73 


0.01 


OGLE SMC 


Max 


0.013 


0.324 ±0.021 


0.013 


0.396 ± 0.039 


0.011 


0.207 ±0.071 


3.16 


0.04 


2.23 


0.03 


Mean 


0.008 


0.272 ±0.017 


0.007 


0.229 ± 0.030 


0.012 


0.340 ± 0.074 


1.77 


0.17 


1.65 


0.10 


Min 


0.007 


0.276 ±0.015 


0.006 


0.282 ±0.027 


0.010 


0.229 ± 0.066 


0.42 


0.66 


0.87 


0.38 



" MRSS denotes mean of residual sum of squares from the fit. MRSSaii, MRSSs and MRSSl are for the fits to all, short and long 
Cepheids in the sample, respectively. 
' Slopes for the period-color relations for all, short and long period Cepheids in the sample. 



the coefficients of the linear regression. Hence the only way 
that inclusion of errors in our least squares fits can affect the 
significance of our F statistic results is if the errors are sys- 
tematically greater for large numbers of either short or long 
period Cepheids. If this is the case, it will affect not only 
our work but also published values of the LMC PL relations 
in the V and I band (e.g., in lUdalski et alJll^9^a^ . 

We now discuss the assumptions required by the F test. 
These are that the error (e) in equation (4) is constant for 
all the stars in our study, and further they are normally 



distributed with zero mean. To test these assumptions, we 
analyze the OGLE LMC data in greater detail. Figure 2] 
shows the plots of residuals from a single PC fit and two PC 
fits to the entire data against log(P). We see no discernible 
trend in the size of the residuals with period and are led 
to the conclusion that the residuals are homoskedastic (i.e. 
a ~ constant). This supports our approach of adopting a 
constant a to gauge photometric errors when doing a least 
squares fit. In addition, the residuals in Figure|3]from the fit 
with a single PC relation show a trend (though not in the 
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Table 2. Amplitude-color relations at maximum, mean and minimum light. 



Phase 






MRSS'^g 


Slope^ 


MRSSl 


Slope\ 


F 


P(F) 


t 


P(t) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


Galactic 


Max 


0.010 


-0.098 ±0.042 


0.007 


-0.313 ± 0.058 


0.012 


-0.169 ±0.088 


15.54 


0.00 


1.50 


0.14 


Mean 


0.010 


0.241 ± 0.042 


0.007 


0.011 ± 0.059 


0.014 


0.275 ±0.095 


10.98 


0.00 


2.65 


0.01 


Min 


0.009 


0.322 ± 0.041 


0.007 


0.145 ± 0.058 


0.013 


0.276 ±0.089 


9.76 


0.00 


1.36 


0.18 


OGLE LMC 


Max 


0.010 


-0.285 ±0.023 


0.007 


-0.393 ± 0.023 


0.007 


-0.160 ±0.058 


93.12 


0.00 


3.62 


0.00 


Mean 


0.011 


0.016 ± 0.025 


0.008 


-0.118 ± 0.024 


0.008 


0.322 ± 0.064 


95.99 


0.00 


6.42 


0.00 


Min 


0.011 


0.140 ±0.025 


0.009 


0.017 ±0.025 


0.009 


0.321 ±0.065 


92.66 


0.00 


4.33 


0.00 


OGLE SMC 


Max 


0.010 


-0.438 ±0.021 


0.007 


-0.419 ± 0.018 


0.010 


-0.308 ± 0.070 


80.45 


0.00 


1.79 


0.07 


Mean 


0.012 


-0.158 ±0.023 


0.007 


-0.144 ±0.019 


0.017 


0.034 ± 0.092 


75.40 


0.00 


2.54 


0.01 


Min 


0.012 


-0.065 ±0.023 


0.008 


-0.047 ± 0.020 


0.011 


0.104 ±0.076 


82.45 


0.00 


2.21 


0.03 



" MRSS denotes mean of residual sum of squares from the fit. MRS Sail, MRSS$ and AIRSSj^ are for the fits to all, short and long 
Cepheids in the sample, respectively. 
^ Slopes for the amplitude-color relations for all, short and long period Cepheids in the sample. 



size of the residuals) in long period Cepheids, which is re- 
duced when using two PC relations for short and long period 
Cepheids to fit the data. In order to test that the residuals 
are normally distributed, we can plot the quantiles of the 
distribution resulting from the ordered residuals against the 
expected quantiles from a normal distribution: a qq-plot. If 
the residuals are indeed normally distributed then this plot 
should be close to the line y = x. Figure |^ shows such a 
plot. There is some small departure from normality at the 
extremes but we contend that this plot justifies our use of 
the F test. 

There is a great deal of information in Table & |51 
but for the context of the present study we summarize the 
results in the following subsections. 

3.1 Results of period-color relations 

For the Galactic Cepheids, the slope of the long period PC 
relation at maximum light is close to zero. Further, the over- 
all slope at maxi mum light is the shallowest. This then sup- 
ports the work of lCodd (Il947f) . In addition, the F and t test 
indicates that the PC relation at maximum light is not con- 
sistent with a single slope. However, the mean and minimum 
light PC relations for the Galactic Cepheids are consistent 
with a single line. 

For the LMC Cepheids, the slopes of the long period PC 
relation at maximum light is also close to zero, decreasing 
significantly from its value for short period Cepheids. The 
slope at mean and minimum light all increase when going 
from short to long period Cepheids. In the LMC, we see 
that the PC relation at maximum light is not consistent 
with a single line at the 93% confidence level. However, in 
terms of the F test, the PC relation at mean and minimum 
light does show evidence to support a break at 10 days. 
Note that at both these phases, the slope of the PC relation 
for longer period Cepheids is much steeper than for shorter 
period Cepheids. 
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Figure 4. Plots of the residuals of LMC period-color relation at 
mean light. Top: (a) The residuals from the long and short period 
PC relations as a function of log(P). Bottom: (b) The residuals 
from using a single PC relation as function of log(P). From the 
plots it is clear that the residuals are homoskedastic. In addition, 
the residuals from the fit with a single PC relation show a trend 
in long period Cepheids. This trend is reduced if two PC relations 
were used to fit the data. 



In the case of the SMC again the PC relation at max- 
imum light is much flatter for longer period Cepheids than 
for shorter period. The relation at mean light for longer pe- 
riods is steeper than at shorter period. Only the PC relation 
at maximum light shows statistical evidence of a break at 
10 days. 

We note that the slope of PC relation at maximum light 
always decreases in going from short to long period Cepheids 
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Figure 5. Plot of the ordered residual quantiles vs. normal quan- 
tiles, the qq-plot. The line is the case oi y = x. 

in all three galaxies. In contrast, the slope at mean and 
minimum hght always increases except for the Galactic and 
SMC Cepheids at minimum light. We also make the inter- 
esting observation that the dispersion of the PC relation, 
whether for the entire sample or either short or long pe- 
riod group, is always the smallest at minimum light in these 
three galaxies. A glance at the F statistic for mean hght for 
the LMC Cepheids shows that we can reject the null hy- 
pothesis of a single PC relation for long and short period 
LMC Cepheids at greater than the 99% confidence level, 
confirming the finding of the broken PC relatio n in LMC 
ifTamm ann fc Reindi2o'o3 : lTammann et alj20o3) . However, 
the Galactic and SMC mean light PC relation is consistent 
with one line. 



3.2 Results of amplitude-color relations 

The AC relations all show statistical evidence of two lines 
at greater than 99.9% confidence level with the F test. At 
mean and minimum light, we see that the slopes of the AC 
relation increase significantly when going from short to long 
period Cepheids in all three galaxies. However at maximum 
light the slope becomes shallower for long period Cepheids. 
Also, the slopes of the AC relation at maximum light are 
always negative for both short and long periods Cepheids in 
Galaxy, LMC and SMC, which we do not see in the mean and 
minimum light. In addition, the slopes are fiat for the short 
period Galactic Cepheids and long period SMC Cepheids at 
mean light, but not for the Cepheids in LMC. 

We also remark on some important differences in the AC 
plane between Galactic and Magellanic Cloud Cepheids. At 
minimum light in the LMC and SMC, the slope of the AC 
relation is 0.017± 0.025 and -0.047± 0.020, respectively, for 
short period Cepheids, whereas it is 0.145 ±0.058 for Galac- 
tic Cepheids - more than 2a difference. For longer period 
Cepheids in the LMC and SMC, this slope becomes positive 
(0.321 ± 0.065 and 0.104 ± 0.076), whilst the longer period 
AC slope in the Galaxy is always significantly above zero. 



Even though Table |5| shows that AC relations are sig- 
nificantly different for short and long period Cepheids in the 
SMC, Figure |3Jb) suggests that the effect is much reduced 
in the SMC as compared to, e.g. Figure 0^b). One possible 
reason for this is that SMC Cep heids have lower amplitude s 
than LMC or Galactic Cepheids dPaczvriski fc Pindoil200oll . 
Equation (1) shows that even if the range of Tmax or Tmm 
is narrow, but if the amplitudes are not large, then Tmin or 
Tmax will not be driven to such low or high values. 

3.3 Combining the PC and AC relations 

As discussed in the Introduction, equation (1) predicts that 
if the PC relation is fiat at maximum light, then there will 
be an AC relation at minimum light, and vice versa. We 
see some evidence to support this idea from Figure [TO and 
from Table Qfc 121 In summary: 

Galactic short period Cepheids: PC relation steep at max- 
imum — > AC ~fiat at minimum. 

Galactic long period Cepheids: PC relation flat at maxi- 
mum AC relations '^steep at minimum. PC relation steep 
at minimum — + AC relation ~fiat at maximum. 

LMC short period Cepheids: PC relation ~steep at max- 
imum — > AC relation fiat at minimum. 

LMC long period Cepheids: PC relation fiat at maximum 
~* AC relations steep at minimum; PC relation steep at 
minimum — > AC relation ~fiat at maximum. 

SMC short period Cepheids: PC relation steep at maxi- 
mum AC relation ~flat at minimum. 



4 CONCLUSION AND DISCUSSION 

We have presented in Figures [ll3l new obse rvational charac- 
teristi cs of Cepheids. Following the work of iTammann et all 
i2002h . we apply rigorous statistical tests to show that at 
mean light, there is a change in the LMC PC relation be- 
tween short (log(P) < 1.0) and long period Cepheids, as 
shown in Figure \^a) . In addition, the LMC data also ex- 
hibit a change in the PC relation at minimum light. How- 
ever, we find no such change at mean and minimum light 
in the Galactic and SMC Cepheids. We also find that the 
PC relations at minimum light generally show the smallest 
scatter among the three galaxies, as compare to the phases 
at maximum or mean light. Following the work of SKM 
we study amplitude-color diagrams at maximum, mean and 
minimum light and find very different behavior in the three 
galaxies for short and long period Cepheids. This difference 
not only occurs within a given galaxy between the short 
and long period Cepheids, but also occurs from galaxy to 
galaxy. For example, the behavior of short and long period 
Cepheids in the LMC is very different (i.e. Figure |5J, and 
the behavior of Galactic and LMC Cepheids of short period 
is also very different. Thus we note that there is an effect 
with both period and metallicity. Further work with state 
of the art pulsation codes is under way to confront these 
observations with model calculations. In addition, the PC 
relations clearly show greater structure than a simple two 
line regression as used here. For example, the LMC data in- 
dicates another break at log(P) ^ 1.2 (see, e.g.. Figure QJ. 
This will be the subject of a future paper. 
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Table 3. F test significance results with color cuts for LMC 
Cephcids. 



(V - I)o 


F 


P(F) 


F 


P(F) 


F 


P{F) 


Range 


Max 


Max 


Mean 


Mean 


Min 


Min 


- 1.9 


2.67 


0.07 


13.2 


0.00 


10.4 


0.00 


- 1.1 


2.95 


0.05 


14.0 


0.00 


10.7 


0.00 


- 1.0 


4.45 


0.01 


10.4 


0.00 


7.86 


0.00 


- 0.9 


3.34 


0.04 


7.69 


0.00 


6.38 


0.00 


0.3 - 1.9 


2.67 


0.07 


13.2 


0.00 


10.4 


0.00 


0.4 - 1.9 


2.48 


0.08 


14.8 


0.00 


12.2 


0.00 


0.5 - 1.9 


2.41 


0.09 


16.3 


0.00 


13.0 


0.00 


0.5 - 0.9 


3.16 


0.04 


10.4 


0.00 


8.76 


0.00 



4.1 Testing of Systematics 

In this subsection we discuss some of the systematic effects 
that may affect our results, as follows: 

(i) Could reddening errors produce the results displayed 
in Figures 11131 '' Consider the bottom panels of Figure 
and IHb) showing AC relations at minimum light for the 
Galaxy and LMC. In contrast to the Galactic counterpart, 
the AC relation for short period LMC Cepheids is fiat. It 
is difficult to imagine a situation where the published red- 
dening for these short period stars are in error to the extent 
of making this relation flat when it should have a slope like 
the Galactic relation. Furthermore, we use the E{B — V) 
values from the literature, which would not only affect our 
results but also other published results that are based on 
these values. 

(ii) Could outliers, perhaps stars that have been misclas- 
sified as Cepheids, produce our significant results? We exam- 
ine the OGLE LMC Cepheid data to investigate this ques- 
tion. The plot of mean [V — I)o color against log(P) in the 
middle panel of Figure |2I a) shows a number of stars which 
seem to have very red colors. We exclude these stars and 
repeat our analysis for PC relations in the LMC. We note 
that it is possible to reduce the significance of our result at 
mean light by excluding stars that are too red. However if 
we exclude stars that are too blue, the F test again becomes 
significant. This conclusion can be seen in Table 01 where 
we apply various color cuts to the LMC data and then cal- 
culate the F test. Excluding the Cepheids with too red or 
blue color does not alter the results we have in Section 3.1. 

(iii) Could we reduce the significance of our results by 
"judiciously" removing certain stars? Since we only consider 
stars with log(P) > 0.4, if we extend this cutoff to 0.6 {P 
= 3.98 days) and repeat our analysis, we get very similar 
results. Also, excluding some of the longer period Cepheids, 
those with periods greater than log(P) > 1.4 increases the 
significance of our results. Hence the period cut we applied 
in this study would not significantly alter the results or con- 
clusions we have. 

While excluding stars based on color or period cuts as 
mentioned above, in some cases, can reduce the significance 
of our F test results at mean light (though never to a P[F) 
value less than 0.1), such experiments have very little effect 
at minimum light. If we accept the assertion that mean light 
properties of Cepheids are an average of Cepheid properties 
at all pulsation phases, then irrespective of mean light. Ta- 



xable 4. F test significance results for the LMC PL relation 



phase 


Fv band 


P(F)v band 


band 


P{F)l band 


Max 


0.32 


0.73 


0.04 


0.96 


Mean 


8.86 


0.00 


6.25 


0.00 


Min 


17.5 


0.00 


19.0 


0.00 



bles I1I2I and Figures I1I3I clearly demonstrate a significant 
change in Cepheid observational properties between short 
and long period Cepheids in the Galaxy, LMC and SMC. 
Further if we accept these exclusions, they would also have 
an effect on the PL relations in both V and I bands and hence 
on the currently accepted extra-galactic distance scale. 



4.2 The period-luminosity relations in LMC 

Since the PL relation is affected by the PC relation, it is 
of interest to examine the PL relation in the OGLE LMC 
Cepheids, which show a break in the PC relation. Do the 
LMC data support a single PL relation or is the data con- 
sistent with a break a t a period of 10 days as shown in 
iTammann et al.l ll2002ll ? We can use the F test to investi- 
gate the LMC PL relations in the V and I band as in Section 
3. The results of the F test are shown in Table 0] where col- 
umn 1 displays the phase and column 2 and 3 show the F 
statistic and its probability, under the null hypothesis of a 
single line, for the V band PL relations. Similar results for 
I band PL relations are listed in column 4 and 5 in Table 
13] From the table, both V and I band PL relations at mean 
and minimum light are significant - that is the data are more 
consistent with a model where the slopes of the V and I band 
PL relation are different for short and long period Cepheids. 
However, the data are consistent with a single slope, in both 

V and I bands, for the PL relation at maximum light, even 
though the corresponding PC relation shows some evidence 
of a break at 10 days (see Table 0. 

The actual slopes of the V and I band PL relation at 
maximum, mean and minimum light for short and long pe- 
riod Cepheids are given in Table |^ The second column in 
this table gives the overall slope. The third and fourth col- 
umn give the short period and long period slopes, respec- 
tively. We clearly see that the slope at maximum light is vir- 
tually identical for short and long period Cepheids in both 

V and I bands. This is very different to the situation at 
mean and minimum light. How can the PC relation at max- 
imum light show evidence of a break at log(P) = 1.0 yet the 
PL relation in V and I be consistent with one PL relation? 
This is certainly not the case for the PL relation at mean 
light in the LMC. This occurs because the LMC PC rela- 
tion at mean light becomes steeper for long period Cepheids 
whereas at maximum light, the PC relation in the LMC be- 
comes shallower for longer period Cepheids. Since the PC 
relation affects the PL relation, this shallow dependence at 
maximum light for longer period Cepheids suggests that the 
PL relation at that phase would not be affected too much. 
A detailed analysis of the PL relations at maximum, mean 
and minimum light will be presented in forthcoming paper. 
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Table 5. Slopes of the V and I band PL relation at maximum, 
mean and minimum light for Cepheids in the LMC. 



Phase 


V slope (all) 


V slope (short) 


V slope (long) 


Max 


-2.74 ±0.05 


-2.73 ± 0.08 


-2.97 ±0.26 


Mean 


-2.75 ±0.04 


-2.95 ± 0.07 


-2.35 ±0.25 


Min 


-2.59 ±0.04 


-2.85 ± 0.06 


-1.84 ±0.26 


Phase 


I slope (all) 


I slope (short) 


I slope (long) 


Max 


-2.95 ±0.03 


-2.96 ± 0.05 


-2.94 ±0.19 


Mean 


-2.99 ±0.03 


-3.10 ±0.04 


-2.94 ±0.18 


Min 


-2.88 ±0.03 


-3.06 ±0.04 


-2.33 ±0.17 



4.3 Future work 

We new briefly summarize other implications from the 
present study. However, the detailed analysis of these topics 
are beyond the scope of this paper, and will be addressed in 
future papers. 

(i) Metallicity dependence on PC relation: From 
the period-mean density relation for a pulsator and the 
Stefan-Boltzmann law, it is easy to show that: log(P) = 
alog(L) + /31og(A/) + 7log(Te//) + const. The parameters 
(a, f}, 7) and the constant term can be determined from stel- 
lar pulsation calculations. For example: 

log(P) = 0.841og(L) - 0.681og(M) - 3.481og(re//) + 
const., from Ivan Albada fc Baker! il973l) . 

log(P) = 0.931og(L) - . 771og (M) - 3.541og(re//) + 
const., from lSimon fc Youn3 1^93). 

log(P) = 0.821og(L) - 0.601og(M) - 3.55 log(re//) + 
const., from ..Beaulieu et al. (.2001). 

In addition, the luminosity and mass should also obey 
an M-L relation from stellar evolution calculations. This M- 
L relation predicts that lo wer metallicity Cep heids will have 
higher luminosity (see, e.g.. lBono et al.l2000h for given mass. 
Hence, at fixed period, lower metallicity Cepheids should 
be hotter tha n higher metallicity Cepheids. For example, 
iLanev fc Stob ie (1986) showed that the long period SMC 
Cepheids are hotter than LMC Cepheids at given period by 
- 200is:. 

Since the Magellanic Clouds (MC) have lower metallicity 
than the Galaxy, C epheids in the MC shou ld be bluer than 
Galactic Cepheids. iTammann et alJ ll2003l) have reported 
this finding in their paper. In this study, we also found that 
the short period MC Cepheids generally have bluer {V — I)o 
color than Galactic Cepheids, for a given period, at maxi- 
mum, mean and minimum light. However, larger errors in 
the PC relations for long period Cepheids make such a state- 
ment more contentious for long period Cepheids. We intend 
to investigate this in detail in a futu re paper. 

(ii) Estimation of color excess : lFerniel il994|) u s ed the 
theory of SKM and the original suggestion of ICodd (Il947l) 
to establish a relation between {B — V) color at maximum 
light, the V band amplitude, period and the color excess 
E{B — V) for Galactic Cepheids. Figure shows that an- 
other approach would be using the properties of Cepheids 
at minimum light to estimate the color excess, since the 
PC relations at minimum light in generally show a smaller 
dispersion (see Section 3.1). Furthermore, the Galactic PC 



relation at maximum light indicates a break at 10 days (Ta- 
ble 0, but not in the PC relation at minimum light. In 
future work, we plan to investigate the possibility of using 
this relation as a way to determine reddening. 
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